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^ ' Abstract 

A structured mixed phase consisting of quark and hadron phases is numericahy 
studied with the Coulomb screening effect and the surface effect. We carefully intro- 
duced the Coulomb potential, so that a geometrical structure becomes mechanically 
! unstable when the surface tension is large. Charge densities are largely rearranged by 

\^ I the screening effect, and thereby the equation of state shows the similar behavior to 

that given by the Maxwell construction. Therefore, although bulk calculations with 
. the Gibbs conditions show that the mixed phase may exist in a wide density region, we 

Q I can see it is restricted to a narrow density region by the surface effect and the Coulomb 

■ screening effect. 

o - 
^. 

1 Introduction 

^ I It has been believed that hadron matter changes to quark matter at high-density region 

by way of the "deconfinement phase transition". Unfortunately the deconfinement phase 
transition have not been well understood up to now, and many authors have studied it by 
model calculations or by first-principle calculations like lattice QCD. These studies are now 
5^ i developing, and many exciting results have been reported. Properties of quark matter have 

been actively studied theoretically in quark-gluon plasma, color superconductivity fP, "5] or 
magnetism and experimentally in relativistic heavy-ion collision (RHIC), HERA or 

early universe and compact stars [HI E] ■ 

When we calculate uniform hadron matter (nuclear matter) and quark matter at zero 
temperature separately, by using the MIT bag model, we can expect the first-order phase 
transition as seen in Figs.QandEl We can see that quark matter is an energetically favorable 
state at high-density region, p > pc (Fig. Q). As we can see in Fig. |21 the thermodynamic 
potential of quark matter becomes lower at higher baryon-number chemical potential. These 
results suggest the deconfinement phase transition at high densities. 

The features of the deconfinement phase transition have not been fully elucidated yet. 
We assume here that it is the first order phase transition, and use the bag model for simplic- 
ity. Then the thermodynamically forbidden region appears in the equation of state (EOS) 
and we can expect the mixed phase, the hadron-quark mixed phase, in some density re- 
gion, which may exist in inner core region of neutron stars and during the hadronization of 
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high-temperature quark-gluon plasma at RHIC experiment. We have to apply the Gibbs con- 
ditions (GC) to get EOS in thermodynamic equilibrium: GC demand chemical equilibrium, 
pressure valance and thermal equilibrium between two phases; 
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where /Xg and /iq are baryon-number and charge chemical potentials, respectively. Thermal 
equilibrium is implicitly achieved at T = 0. Note that there are two independent chemical 
potentials in this phase transition. In such a case the system should be much different from 
the liquid-vapor phase transition which is described by only one chemical potential. 
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Figure 1: Energy density e for uniform 
hadron matter and quark matter. Uniform 
quark matter is energetically favorable in 
high-density region, while hadron matter in 
low-density region. 



Figure 2: Thermodynamic potential of 
uniform hadron matter and quark matter as 
a function of baryon chemical potential fi-Q- 
Uniform quark matter is energetically favor- 
able in high /iB region while hadron matter 
in low /iB region. 



These GC must be fulfilled in the hadron-quark mixed phase. On the other hand, the 
Maxwell construction (MC) may be very familiar and has been used by many authors to get 
EOS for the first order phase transitions [HI El EH E] • It is well known that MC is a correct 
prescription to derive EOS for the liquid-vapor phase transition. However, Glendenning jj^ 
pointed out that MC is not appropriate for the hadron-quark mixed phase: one of GC about 
the charge chemical equilibrium, /ig"^'^'^ = /ig^'^™'^, is not satisfied in MC, since the local 
charge neutrality is implicitly assumed without imposing this condition. He emphasized 
that the local charge neutrality is too restrictive and each hadron or quark phase may have 
a net charge because only the total charge must be kept neutral. 

When we use GC in the bulk calculation, which we explain in detail later, we can see 
the mixed phase appears in a large density region. The pressure is not constant as density 
changes in the mixed phase, while only a constant pressure is obtained from MC as shown 
in Refs. 

However, the bulk calculation is too simple and Heiselberg et al. ^3] claimed the im- 
portance of including the finite-size effects, i.e., the surface tension at the hadron-quark 
boundary and the Coulomb interaction energy. They studied the quark droplet immersed in 
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hadron matter and found that it is energetically unfavorable if the surface tension is large 
enough. On the contrary, however, if the surface tension is not large the mixed phase can 
exist in some density region. 

Glendenning and Pei suggested the crystalline structure by a bulk calculation using 
the small surface tension: one phase is immersed in another phase with various geometrical 
structures; "droplet" , "rod" , "slab" , "tube" , and "bubble" . These are called the structured 
mixed phases (SMP). Applying the results based on the bulk calculation to neutron stars, 
they suggested that there could develop SMP in the core region for several kilo meters in 
thickness. 

At a first glance, 
this view seems to be 
reasonable and there 
may appear the mixed 
phase in a large den- 
sity region. However, 
there are still many 
points to be elucidated 
about the finite-size ef- 
fects. Voskresensky et 
al. emphasized that the 
proper treatment of the 
Coulomb interaction is 
important in the mixed 
phase |16J. 

First note that there 
is the relation between 
chemical potential and 

the Coulomb potential by way of the gauge transformation: chemical potential is not well 
defined before the gauge fixing. Secondly, the charge chemical equilibrium could be rather 
satisfied even in MC, once the Coulomb potential Vcoui is incorporated: 
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Figure 3: Schematic view of the structured mixed phases. 
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whereas 

quark r^quark / hadron T/hadron fo\ 
^ '^Coul T f^Q ~ t^Coul ■ l-JJ 

Note that the electron chemical potential fig = fiq and the electron number density can 
be expressed in terms of the combination, (/iq — Vcoui) (see Eq. (48)), when the Coulomb 
potential Vcoui is introduced. Therefore, Eqs. © niean that the charge chemical poten- 
tial is equal, while the electron number is different between the hadron and quark phases. 
Equation Q is reduced to fif^^^^ ^ in the absence of the Coulomb potential, which is 

the previous claim that the charge chemical equilibrium is apparently violated in MC. Thus 
we can see that the previous claim means nothing but the difference of the electron number 
in two phases fH]- Thirdly, it is important to take into account the Coulomb screening 
effect^. We have to solve the Poisson equation consistently with other equations of motion 



^The Coulomb screening effect on kaon condensation is studied by Norsen and Reddy [T7). 
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for charged particles. Actually, Voskresensky et al. showed SMP is mechanically unstable by 
the Coulomb screening effect when the surface tension is not small. If SMP is mechanically 
unstable, the phase transition should be similar to that in MC. 

It can be easily seen that there should occur rearrangement of charged particle densities 
by the Coulomb interaction. On the contrary such rearrangement modifies the Coulomb 
potential. As a result the Poisson equation becomes highly non-linear and difficult to solve 
analytically. In their study, a linear approximation (RPA) was employed to solve the Poisson 
equation analytically. It is valid to use the approximation for pointing out the important 
property of SMP, but it may be conceivable that various charge properties like "global charge 
neutrality" for GC, "local charge neutrality" for MC, and "the Coulomb screening effect" 
have important roles in the mixed phase. Therefore, solving the Poisson equation without 
any approximation is of much significance. We have reported a preliminary result for the case 
of droplet [1B|. It would be also very interesting to derive EOS in our consistent calculation 
for the hadron-quark matter phase transition. 



2 Bulk calculation with the Gibbs conditions 

Now we consider two infinite matters separated by a sharp boundary: uniform quark 
matter and hadron matter. We consider the quark phase consists of u, d, s quarks and 
electron, and the hadron phase proton, neutron and electron. We discard the Coulomb 
interaction in the calculation. Then we can evaluate the total thermodynamic potential 

f^tot = + fid + + fin + fip + fie- (4) 

The explicit expressions of fij for non-uniform matter are given in the next section. We use 
them by replacing the space dependent quantities by constants for each uniform matter. 

Introducing the volume fraction of the quark phase /, we impose the global charge 
neutrality: the total charge density vanishes, 

/pQ + (1 - /)pH = 0, (5) 

where pq and pu are the net charge densities of the quark and hadron phases, respectively, 

PQ = - ^Pd - - Pe, 

PH = Pp - Pe- 

Note that we have now six chemical potentials; p^, pd, Ps, Pp, Pn{^ A^b), A*e- We first 
consider (3 equilibrium in each phase and chemical equilibrium at the hadron-quark boundary. 
Thus all chemical potentials are constant. The chemical equilibrium conditions then are: 

Pu + Pe= Pd, (7) 

Pd = Ps, (8) 

in the quark phase, 

Pp + Pe = Pn (9) 
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in the hadron phase, and 



/^n = /Uu + 2/id, (10) 

/ip = 2/iu + /id, (11) 

at the hadron-quark interface. The last condition (lllj) can be derived from other four condi- 
tions, so that there are left four independent conditions for chemical equilibrium. Therefore, 
if we give two chemical potentials /iB and fi^ by hand, we can determine these four chemical 
potentials; fi^, /id, /is and /Xp from these four equations. 

Next, we can determine /ie by the global charge neutrality condition ©• / is still unknown 
at this point and finally, we find the optimal value of / by using one of GC; p^^^^^ = phadron^ 
where pressure P is given by the thermodynamic relation: P = —Q/V. Thus once /ie is 
given, all other values /ij (i = u, d, s,p, e) and / can be obtained. 

We show the results of the bulk calculation. We present the charge densities in Fig. EJ 
We can see the total charge density is zero, while the quark phase is negatively charged and 
the hadron phase positively charged. Note that MC always gives a null charge density in 
each phase due to the local charge neutrality [TSj. 

Figure El shows pressure versus baryon- 
number density. The most important differ- 
ence between "Maxwell" and "Bulk Gibbs" 
is that the pressure given by MC is constant 
while that given by the bulk calculation with 
GC is density dependent. Remember that the 
outstanding feature in these two results comes 
from that each hadron or quark phase has a 
net charge in "Bulk Gibbs" , while each phase 
is neutral in "Maxwell" . 

The phase diagram in the /xb — /^e plane is 
presented in Fig. IHl We can easily see that the 
charge chemical potential /ie is much different 
between uniform hadron matter and quark 
matter for MC, which means chemical equi- 
librium condition for is not satisfied in MC. Thus it could be said that MC is not a proper 
way in the bulk calculation. Glendenning J2] pointed out this defect in MC and showed 
another result (Bulk Gibbs) within the bulk calculation. One can easily see that the mixed 
phase appears in a wide /ab region in Fig. El 

After the claim by Glendenning, Heiselberg et al. demonstrated that the finite- 
size effects may disfavor the mixed phase by extending the bulk calculation to include the 
Coulomb interaction and the surface energy. On the other hand, Glendenning and Pei jT^j 
suggested "crystalline structures of the mixed phase" because SMP are energetically more 
favorable than uniform matter even including the surface and Coulomb effects. 

Voskresensky et al. pointed out the importance of the Coulomb screening effect and 
suggested the usefulness of MC for the mixed phase. We have seen MC is apparently incorrect 
because if we apply MC to this mixed phase, one of GC, charge chemical equilibrium, is 
violated. However, once the Coulomb interaction is taken into account, there is still a 
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Figure 4: Charge density of each phase. The 
total charge density is always zero. 
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Figure 5: Pressure of uniform matter, that 
given by the bulk calculation with GC (Bulk 
Gibbs) and that given by MC (Maxwell). 



Figure 6: Phase diagram in the chemical 
potential plane. H means uniform hadron 
matter and Q uniform quark matter. Other 
notations have the same meaning as in Fig. 



room to satisfy the condition without spoiling the picture given by MC, as noted in the 
previous section. Since the Coulomb screening effect had been taken into account with an 
approximation in their study, we may still address a question, "is the picture derived from 
MC is completely meaningless when we include the Coulomb interaction in a proper way?" . 
All we have to do next is to solve the Poisson equation without any approximation, strictly 
satisfying GC. 



3 Formalism 



3.1 Density Functional Theory 

We use the idea of the Density Functional Theory (DFT) jTHl ED] to study the quark- 
hadron mixed phase. We employ here the Thomas-Fermi approximation, which gives the 
simple expressions for the energies. We summarize below some energy expressions with the 
local density approximation (LDA); the expressions are first derived for uniform system, and 
then we regard densities as the space-dependent functions. 

First, the kinetic energy is simply expressed as 



e(p(r)) 
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where m is the particle mass and the Fermi momentum Pf(^) is expressed by the density 

profile p{r) as ppir) = {'n'^p{r)Y . 

Next, we consider the interaction energy of quarks. We take the first-order contribution 
which comes from the one-gluon exchange interaction in uniform quark matter, 



Here Aq is the SU(3) Gell-Mann matrix and D^^ the gluon propagator [21]. By way of the 
Wick contraction, the Hartree and Fock terms are derived: 









(14) 



for the Hartree term, and 



^(g)(-^^77/.y )V^(^)^'^^(g-^)^(A:) 



q 
k 

q 



/\ a q-k 



(15) 



for the Fock term. Due to the traceless property of A^, TrAa = 0, the Hartree term gives a 
null contribution in the color-singlet quark matter. 



_l 7 

%l){k) {-igiu^ Hik) =Tr 



0. 



(16) 



Therefore only the Fock term contributes to the interaction energy, which becomes 
k 




Tf^T^^Tr {l,Gn{q)l.Gnm D^^{q - k) (17) 
(2vr) (27r) 



£q d^k 1 2 . \ 



2 — ^Fock- 



(27r)3(27r)3 ' " Ak\^ ^-"^ " [k - qf 
Here Gd is the density- dependent quark propagator The subscript / denotes the flavor. 



mj the quark mass, k^ = [\k\" + rrij^ 



1 /2 

uo _ (\7^\2 I g^j-^^ _ ()(^kpj, — j/gj^ with the Fermi momentum 

kpf. After some manipulation we easily find 
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(19) 



where xj = ^^^^ ^ and rjf = . /l + x'i, and we use the relation PF{{r) = (7r^pj(r)) . 
nij V 

Equation (fT^ becomes epock = iff S/Pf/ massless particles. 

For the interaction energy of nucleons, we use for simplicity the effective potential 

parametrized by densities [TB] . 

Poir) 

2 



(Pn(y) +Pp(r)) / pn(r) + Pp _ 

+aat (Pn(r) + pp(r)) |^ /^"(^)+/^p(^) _ 1^ , (20) 

where, 5*0, Kq, ebind, and Csat are adjustable parameters to satisfy the saturation properties 
of nuclear matter. 



3.2 Thermodynamic potential 
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nucleon quark 

Figure 7: Wigner-Seitz ap- 
proximation. 



We consider a system which consists of the hadron and 
quark phases. We divide the whole space into equivalent and 
charge- neutral Wigner-Seitz cells with the size Rw and the 
droplet size R as illustrated in Fig. [7| Here we only consider 
the droplet phase. Following DFT we begin with the ther- 
modynamic potential in terms of particle density profiles and 
chemical potentials, p^ and pf , 

^tot = E{p,{r)) -Y.fi? drpfir) - $^pf / drpf{r). 

i Jo i Jr 

(21) 

To take into account the Coulomb interaction, we write it in 
terms of the i-particle density profiles; pi{r) with Qi being the 
particle charge (Q = — e < for the electron). 



Applying the Laplacian on Eq. (j^ . we can easily de- 
rive the Poisson equation V^Vcoui = 47re^ Qjpj(r). The 
Coulomb interaction energy Ey is then expressed as 



V 



(23) 



Including the surface term, the total energy is expressed as 

rR pRw r 

E{pi{r))= / dreQ + / rfren + / dSes + Ey, 
Jo Jr Js 



(24) 
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where es means the surface-energy density, S boundary area, and eg and en the energy 
densities in the quark and hadron phases. Since we poorly know the details of the hadron 

quark interface, we simply approximate the surface energy as / dSes = crS by using the 



surface tension a. Each chemical potential is derived by the equation of motion - — ^ = 0, 

P^ = ^-f^-N,Virl iV, = ^, (25) 
dpi{r) e 

rR rRw 

-E'kin+str = / dreQ + / dren- (26) 
Jo Jr 

As we have seen in the previous section, each phase may have a finite net charge. In 
the previous studies of SMP, the Coulomb interaction between charged particles has been 
treated rather simply: the Coulomb energy was added to the total energy by using the 
volume fraction / and constant densities derived for two infinite matters P21I221- Thus the 
Coulomb screening effect and rearrangement of charge densities are completely discarded. 
One important point we would like to address here is, when there is the Coulomb interaction, 
the gauge variance of chemical potentials should be taken into account. Differentiating the 
expression of chemical potential (j25|l with respect to the Coulomb potential V{r), we get 
the relation [THj, 

= (27) 
AijBjk = Sik, (28) 

where matrices A and B are defined as 



6pi6pj ' dpj 



A,^ B,,^^. (29) 



From these equations, the gauge invariance relation is derived. 

We can immediately see that the chemical potential is gauge variant from this equation: 
when the Coulomb potential is shifted by a constant value, V{r) =^ V{r) — VQ, the chemical 
potential should also be shifted pi =^ pi + A^jVo to keep the gauge invariance. Note that 
the charge density should be observable and gauge invariant; we shall see the constant shift 
of the Coulomb potential is compensated by the redefinition of the chemical potential in 
the expression of the charge density. The Poisson equation is also gauge invariant because 
it is not changed by the constant shift of the Coulomb potential. Therefore, if the Poisson 
equation is not properly taken into account, the gauge invariance is obviously violated. 

Next, we present the explicit form of the thermodynamic potential for each particle. First, 
we consider the quark phase which consists of u, d and s quarks and electron. According to 
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the small current quark mass, we treat u and d quarks as massless particles and only s quark 
massive (150 MeV) in this phase. Electron is also treated as a massless particle. Next, we 
consider hadron phase which consists of proton, neutron and electron. Nucleons are treated 
as non-relativistic particles. The total thermodynamic potential becomes 

Otot = fig + f^H + Es. (31) 

Here Qq and Q-^ are the thermodynamic potentials in the quark and hadron phases, and Es 
the surface energy. 



3.2.1 Quark phase (u, d, s quarks) 

The thermodynamic potentials in the quark phase are given by Eqs. fll2|) . ()19|). (jUJ): 



fiu + ild + + fi!?m + / d'^B (fi : bag constant). 
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es(Ps(?')) - /^sPs(r') + -^Coui(r')ps(r') 



(VV (r\? + 



l^epeir) + Vcoui{r)peir) 



(32) 
(33) 

(34) 
(35) 
■ (36) 



Here, summarizes the electron and another Coulomb contributions. The energy density 
of s quark es{ps{r)) is a function of Ps{r), which is explicitly expressed by using Eqs. (fT^ 
and (fT^ as 



es(Ps(r)) 



3mf 



PFs[r) 



PFs[r} 



Wis 



PFs{r) 
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(37) 
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3.2.2 Hadron phase (non-relativistic nucleons) 

The thermodynamic potentials in the hadron phase read, 

Oh = On + Op + 0^^+ / drepotir) here epot(r) is Eq. (EDI)- 

Jr 



Oh 


= On + 






On 






Jr 




rRw 


Op 
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A 2. 

(Svr^) 3 pi (r) - /in (pp(r), pn(r')) Pn(r') 



(Stt^) >^ (r) - pp (pp(r),pn(r))pp(r) - Vcoui(r)pp(i 



we 



47r2 



(38) 
(39) 
(40) 

.(41) 



Here O^ summarizes the electron and another Coulomb contributions. 

3.3 Equations of motion 

We get the expression of chemical potentials and the Poisson equation from the equation 
of motion ^ = 0, where 0^ = p„(r), pd(r), p,(r), pp(r), p„(r), pe(r), Vcoui(r*)- The Poisson 
equation is explicitly written as 



^Pu[r^ 3 



lpd{r) - ^p,(r)^ 9{R - r) + Pp(r)^(r - R) - pe(r) 



(42) 



The chemical potentials for quarks in the quark phase are derived as 

Pn = 



/id 

Ps 



l + ^)vrU(r)-fWr'), 
2a \ 2 i 1 



eps + 



2a^ 

37T 



m: 



PFs(r')-3^1n 

eps 



eps +PPs(^'j 



m. 



+ gK^ouKr). 



(43) 
(44) 
(45) 



Here, we use eps = \/ ml + p^s ^^"^ Pfs = (7r^Ps(?'))3 . The nucleon chemical potentials in the 
hadron phase and the electron chemical potential are 



Pn 



pL , 2^0 (Pn(r-) - Pp(r)) , f Pnjr) + Ppjr) 

1 H Cbind H — I J- 



2m 



Po 
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Po 



^Ko / pn(r) +pp(r) _ A _^ 2C /"(^) +/^p(^) _ c 
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~ 2m + 2m 



pI^ , pIp 45o(pB(r)-2pp(r))' 



Po 



Vcouiir) 



Pe 



(37r2pe(r))^ + V^Coui(r). 



(46) 
(47) 

(48) 
(49) 
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Note that, e.g., the electron density profile Pc{r) is expressed as 



_ ip.-VcUr)) ^ 

in a gauge invariant fashion (cf. Eq. Q). We impose the (3 equilibrium in the quark and 
hadron phases and the chemical equilibrium at the quark-hadron boundary as in Eqs. ((Tj)- 
flll|) . The pressure contribution from the surface tension is derived as = (^-^ with Vq 
being the volume of the quark droplet. From Eqs. (jH^ and (IHHj). the pressure of each phase 
is expressed as 

P^ = -^, P« = -^, (51) 

with Va being the volume of hadron phase. Therefore the pressure balance condition becomes 

pQ = pH ^ P^, (52) 

which gives the droplet size R. Finally the cell size Rw is determined by the minimum 
condition for fitot- Particle density profiles and R, Rw are completely determined for given 
fiB- Thus we can solve these coupled equations of motion consistently with GC. Note that the 
Coulomb potential is included in a proper way and Vcoui(^) appears in almost all chemical 
potentials. The Coulomb potential is a functional of the charged-particle density profiles and 
in turn densities are functions of the Coulomb potential. As a result, the Poisson equation 
becomes highly non-linear. Since it is difficult to solve analytically, we solve it numerically 
without any approximation, whereas the linear approximation has been used in the previous 
work [TT)] . 



4 Numerical results 

We show a case of quark droplet embedded in hadron matter for the volume fraction 
/ = (i?/i?w)^ = 1/100 and the surface tension cr = 60 MeV/fm^. We can see the Coulomb 
screening effect in Fig. |H1 There is always a minimum in the thermodynamic potential with 
respect to the droplet radius R without the Coulomb screening, due to the balance between 
the Coulomb energy and the surface energy. However, once the Coulomb screening is taken 
into account, the minimum disappears, which shows the mechanical instability caused by the 
Coulomb screening effect. If the surface tension is large enough, the geometrical structure 
(droplet) becomes mechanically unstable by the Coulomb screening. 

Next, we show the density profiles. We can see the fiat density profile of particles in the 
left panel of Fig. IHl where there is no Coulomb screening effect. In the right panel, on the 
other hand, we can see the rearrangement of the charged particles as the Coulomb screening 
effect. 

Let us consider the features of these density profiles in detail. In the quark phase, the 
densities of the negatively charged d, s quarks and electron are reduced, while that of the 
positively charged u quarks are enhanced by the Coulomb screening effect. Remember that 
the quark phase has a negative charge and the hadron phase a positive charge in the hadron- 
quark mixed phase. 
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For the hadron phase, we can see the oppo- 
site effect by the Coulomb screening: the elec- 
tron number increases and the proton number 
is reduced. We see that the charged particles 
in one phase are also affected by the charge 
of another phase: negatively charged parti- 
cles in the quark phase are attracted to the 
surface, while positively charged particles (u 
quarks) are repelled from the surface. On the 
contrary, the positively charged particles in 
the hadron phase (protons) are attracted to 
the surface, while electrons are repelled. Al- 
though the rearrangement effect is not enough 
to establish the local charge neutrality, we can 
understand that the Coulomb screening effect 
is to reduce the local charge in each phase. 

To draw the phase diagram of the hadron- 
quark mixed phase, we have to minimize the 
thermodynamic potential with respect to the 

volume fraction, i.e., we have to change the droplet radius R and Wigner-Seitz cell size R 
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Figure 8: Radial dependence of the thermo- 
dynamic potential Vt for the droplet struc- 
ture. Solid and dashed hues are f2 with and 
without screening, respectively. 
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Figure 9: Density profiles pi without the Coulomb screening effect (left panel) and with the 
Coulomb screening effect (right panel). This is the case of p-Q = 1232 MeV and the volume 
fraction / = {R/ Ry^)^ = 0.01. Each density uniformly spreads in each phase and Vcoui is 
constant ( =0 ) in the left panel, while there is rearrangement of the charge densities and 
Vcoui is spatially dependent in the right panel. 



We show the total thermodynamic potential for the droplet phase in Fig.^J Comparing 
with uniform hadron matter and quark matter, the droplet phase takes the smallest value 
of thermodynamic potential between a certain range of the baryon-number chemical poten- 
tial. To clarify the difference between uniform matter and the droplet phase, we show the 
difference of thermodynamic potential between them in the right panel of Fig. ^| 

We can see a large difference between our result (cr = 40 MeV/fm^) and that given by 
the bulk calculation in Fig. ^| 
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Figure 10: Thermodynamic potential Q (left panel) and its difference from uniform hadron 
matter and quark matter SQ (right panel). The dashed line "Bulk Gibbs", in the right panel, 
shows the difference between the bulk calculation and uniform matter. Solid lines are the 
cases for the droplet phase with a = 10 and 40 MeV/fm^. 



If the difference 6Q is negative, the mixed phase is an energetically favorable state. On 
the contrary, if 6fl is positive, the mixed phase is an unfavorable state. In the right panel 
of Fig. Uni the point given by the Maxwell construction is specified by a circle (/iB ~ 1257 
MeV, 5fi = 0), where the following relations are maintained. 



quark 



hadron/ 



pquark _ p 



hadron 



yquark rp 



hadron 



(53) 



Note that if we naively apply MC, we have to discard one of GC, fi"^^^^^ = fi^^'^'^°'^ in the 
absence of the Coulomb potential. The curve of "Bulk Gibbs" shows that the mixed phase is 
energetically favorable in the wide /iB region, as is already seen in Sec. 2. The difference 6fl 
for "Bulk Gibbs" is up to 10 MeV/fm^, while for the droplet phase {a = 40 MeV/fm^), it is 
only less than 2 MeV/fm^. This shows that if we use the larger value of the surface tension, 
the mixed phase gets more unfavorable. This feature is similar to that reported by Heiselberg 
[TIj or Alford We can see the region of the energetically favorable mixed phase in the 
right panel of Fig. the region becomes narrower than "Bulk Gibbs" , which means that 
the property of the mixed phase is closer to that given by MC. Thus we have seen that the 
mixed phase is very much affected by the Coulomb screening effect and the surface effect, by 
strictly keeping GC. For the larger value of the surface tension, MC is effectively useful in 
the description of the hadron-quark mixed phase. Note that we never violate the condition 
for the charge chemical equilibrium /ig'^'^''°° = /ij^"^'''^ to get these conclusions. Remember 
that the meaning of the condition in the presence of the Coulomb potential is different from 
that in the absence of the Coulomb potential. 



5 Summary and concluding remarks 

In this note we have first seen how the bulk calculation is performed for the hadron- 
quark mixed phase by applying the Gibbs conditions to the system consisting of two infinite 
matters. It gives a wide density region for the mixed phase. Based on this result, some 
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authors also suggested the structured mixed phase with various geometrical structures, by 
including the finite-size effects, the surface and Coulomb energies |22] ■ 

We have examined the finite-size effects in the mixed phase by numerically solving the 
equations of motion for the particle densities and the Poisson equation for the Coulomb 
potential. Our framework is based on the idea of the density functional theory, which, 
we believe, is one of the best theories to treat the structured mixed phases. We have 
demonstrated, by taking the droplet phase as an example, that the Coulomb screening 
effect and rearrangement of the charge densities play an important role for the mechanical 
instability as well as the energy of the mixed phase. As a result we have seen that the 
region of the mixed phase is highly restricted by the Coulomb screening effect as well as 
the surface energy. We have also seen that EOS gives the similar behavior to that given 
by the Maxwell construction, whereas the Maxwell construction is apparently incorrect in a 
system with more than one chemical potential: the Maxwell construction is effectively useful 
in the description of the mixed phase, even in this case. As another case of more than one 
chemical potential, kaon condensation has been also studied j23| and the result is similar to 
the present case. 

We have included the surface tension, but its definite value is not clear and many authors 
treated it as a free parameter. There are also many estimations for the surface tension at 
the hadron-quark interface in lattice QCD [^^3, in shell-model calculations [23 123 I2H| 
and in model calculations based on the Dual-Ginzburg Landau theory j^Hl- If we have the 
realistic value of the surface tension, we can reasonably bring out SMP in the hadron-quark 
phase transition. As these system corresponds to neutron star matter, we have seen that 
the mixed phase should be narrow by the finite-size effects. Our result would restrict the 
allowed SMP region in neutron stars which is suggested by Glendenning ^] • It could 
be said that they should change the property of neutron stars especially in equation of state 

We have assumed in relation to phenomenological implications here that temperature 
is zero. It would be much interesting to include the finite-temperature effect. Then it is 
possible to draw the phase diagram in the /ib — T plane and we can study the properties of 
the deconfinement phase transition. In this study we have used a simple model for quark 
matter to figure out the finite-size effects in the structured mixed phase. However, it has 
been suggested that the color superconductivity is a ground state of quark matter [TJ ■ 
To get more realistic picture of the hadron-quark phase transition, we will need to take into 
account color superconductivity. In the recent studies the mixed phase has been also studied 
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